1 



Decoherence of a Spin Qubit Coupled with a Spin Chain 

Toshifumi Itakura *) 
Department of Physics, Kyoto University, Kyoto 606-8502, Japan 

(Received October 25, 2004) 

In this report, we examine the relaxation phenomena of a spin qubit coupled to a spin 
chain with a 1/r 2 interaction. 

§1. Introduction 

Among various proposals for quantum computations, quantum bits (qubits) in 
solid state materials, such as superconducting Josephson junctions, 1 ) and quantum 
dots, 2 ) -4 ) have the advantage of scalability. Such coherent two level systems consti- 
tute qubits and a quantum computation can be carried out as a unitary operation 
appied to multiple qubit systems. It is essential that this quantum coherencre be 
maintained during computation. However, dephasing is difficult to avoid, due to the 
system's interaction with the environment. The decay of off-diagonal elements of 
the qubit density matrix signals the occurrence of dephasing. Various environments 
can cause dephasing. In solid state systems, the effect of phonons is ubiquitous. 5 ) 
The effect of electromagnetic fluctuation has been extensively studied for Josephson 
junction charge qubits. 6 ) 

The fluctuations of the nuclear spins of impurities can also be a cause of dephas- 
ing. It has recently been shown experimentally that the coupling between the spin of 
an electron in a quantum dot and the environment is very weak. 7 )' 8 )' 10 )' n ) For this 
reason, the dephasing time of a spin qubit conjectured to be very long. However, 
both donor impurities and nuclear spins in semiconductors 9 ) have been suggested as 
possible building blocks for feasible quantum dots architectures. The proposal based 
on experimental findings of quantum computation using an Si 29 array is another 
possibility 9 ^ For a spin qubit system in 2D, the coupling between nearest neighbor 
spin chain can cause a dephasing. When an interaction is between the qubits them- 
selves, it can in principle be incorporated into the quantum computer Hamiltonian, 
although this would lead to more complicated gate sequences. Therefore it is in- 
structive to analyze the error introduced by ignoring some of these interactions, as 
done in the case of dipolar coupled spin qubits. 12 ) For a quantum spin chain with 
a 1/r 2 interaction, an exact expression of the dynamical correlation function has 
been obtained. 14 ) Using this expression, we consider the case in which one qubit is 
coupled to each spin of the spin chain. We examine the relaxation phenomena of 
a spin qubit in array that is coupled to a spin chain with long-range interactions. 
Integrating over the spin chain variables, we obtain the influence function of the 
qubit system. Using this influence function, we examine the dephasing of the qubit's 
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density matrix. In the present study, we especially concentrated on the correlation 
effect between the spin of the qubit and that of the spin chain, and show how the 
characteristic nature of environment appears. The result is that the dephsing rate 
increases as a function of the intrachain interactions, because the fluctuations of the 
spinon is suppressesd interaction increases. 

§2. Hamiltonian 

We examine the Hamiltonian which is given by, 

7~t = Hqb + H qb _ spin + 77 spin (1) 

H qb = HAI Z , (2) 
flint = lNh 2 (A zz I z S z0 + A ± (I + S- + IS+o)), (3) 

N 
i,j=N 

where N is the number of sites, and d(n) = (N/ir) sm^nn/N) — > n as N — > oo. Also 
i is the index of N lattice. A is the effective magnetic field applied to a spin qubit, 
and A zz jnH is the magnitude of the fluctuations in the coupling constant. A system 
of this kind has also been examined in the context of a giant spin that is coupled to 
a surrounding nuclear spin interacting bilinearly with a spin bath in a magnetic field 
and with internal interactions of canonical form. 15 )' 16 ) In this study we examine the 
dephasing behavior by using an exact dynamical correlation function and examine 
the characteristic nature of a bath that interacts with itself. 

§3. Influence function 

We examine first a pure dephasing event (A± = 0), and then we study dephasing 
with a dissipation event (A j_ / 0) . The density matrix of spin qubits is given by 

p(I f z+ ,I f z _) = / [d/ z+ ][d/ z _]exp(-(J 9b [7 z _]-/, b [/, + ]))F[/ z+ ,7 z _], 

(5) 

where I q b[I±] = J hAI z ±, and the influence function (IF) is defined as 

F[I Z+ ,I Z _] = j [dS +i ][dS_ i ]<5(S +i (t)-S_ i (*))p(S i (0),S j (0))exp{l(7[5 + ]-7[5_])}^ 

(6) 

For the following investigation, we assume A±=§ and that the environment system 
action is given by 7[S] = 7d[S] + 7j n t[7 2 , S z ], with free part 

Io[S] = {Azz 2f ? jT l^dhdh^Siih^ih^^ih) 

.6 f , f* , ,dS dS, 
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where 6 = 2nn and Aoo p (ti, i, t2, j) is the free propagator of the environmental 
system at zero temperature, which is defined on a closed time path and has four 
components: 17 ^ 



*00p 



- if a t \ - ( A++(i,ti,i,t 2 ) A+ (i,t u j,t 2 ) \ 

(8) 



When A± = 0, Aoo P — ► ^oop> OIU y z component of the qubit is important. 
For the Berry phase term, because n is odd, the environment is a one-dimensional 
spin-1/2 system. 18 ^ In the present model, the spin chain is a gapless solvable model 
(a special point in spin chain systems). Therefore, we can completely integrate 
over the spin chain degree of freedom and the Berry phase term does not make 
system to be distempered. Thus, introducing the incoming interaction picture for 
the environment, and integrating out the spin chain environment, we can easily verify 
that Eq. © becomes, 17 ) 

F[I Z+ , I z -\ = exp[— i ^2 J J dtidt 2 

(/, + (t 1 )Z\^+ + (0,ti,0,t 2 )/, + (t 2 )+/ 2 _(ti)Z\g -(0,t 1 ,0,t 2 )/^(t 2 ) 
- 4 + (i x )^+-(0, ti, 0, t 2 )I z _{t 2 ) - I z -(h)A Z w + (0, ti, 0, t 2 )I z+ (t 2 ))]. 

(9) 

For convenience we change the coordinates, rj = [I z + + = {h+ — Iz-)/%- 

The coordinates of r\ and £, given by r\ and £ are called the "sojourn" and "blip". In 
terms of these variables, the density matrix is described by p{j] = 1) = | |)(| \,p(rj = 

-i) = 1 1)U U£ = i) = I T)U \,p(Z = -i) = I IXt I- 

Then, the IF becomes 

rti 



F[ V ,t]=eM-i iAzZ T h) [ [ 1 dttdh^G^tt, 0,t 2 ,0)r/(t 2 ) 
4 Jo Jo 

+ V (t 1 )G A (t 1 ,0, t 2 , 0)£(t 2 ) - £(*i)G K (*i, 0, i 2 , 0)£(t 2 )}], (10) 



where G R (t 1 , 0, t 2 , 0), G^ti, 0, t 2 , 0) and the G^(ti, 0, t 2 , 0) are the retarded Green 
function, advanced Green function and Keldysh Green function for i = j = 0. 
It should be noted that the starting Hamiltonians are those given Eqs.([T]) ~ dH). 17 ) 
Because the propagator of a qubit has no time dependence, the blip state and sojourn 
state do not change in time. Therefore, when we choose the initial condition of the 
qubit density matrix to be a coherent state, such as p(t = 0) = ±(| j)(J, | ± | J.)(f |), 
time evolution occurs only in the off-diagonal channel. In addition, even if we start 
from an off-diagonal state, the interaction Hamiltonian does not allow a spin flip 
process. Hence the blip state does not change in time. Therefore we can set £(t) = 
£(= ±1) and rj(t) = for all t. This leads to the following exact expression of the 
dephasing rate : 

p(£,t) = exp[iAt - i ^ AzzlNh ^ f f 1 ' dt 1 dt 2 G K (t 1 ,0,t2,0)] 

4 Jo Jo 
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= exp[iAt _ ( A ^Nh) 2 r l n dhdt2{S z {ti) ^ s'fo)}]. (ii) 

4 Jo Jo 

Here, {Sq,Sq} is the symmetrized correlation function, defined by {A,B} = AB + 
BA, because the Keldysh Green function is a symmetrized correlation function for 
the base of the present Hilbert space. 18 ) 

§4. Results 

First, we examine a pure dephasing event, (i.e., A±=0) Then using analytical 
expressions, the time evolution of the off-diagonal density matrix is given by 

- lnp(£ = ±1, t ) = iAt + {Azzl jf )2 £ jf* dhdhiS^S^h)}. 

Here, E = £(Ai.A 2 ) = ^(Af + \\ - 2AfA2), where v is the spinon velocity. We 
carried out the numerical estimation at zero temperature. 

In the initial regime, ^ < 1, we have -Relnp(t) = i^±M^f-t 2 . This is Gaus- 
sian behavior which is displayed in the initial regime. This comes the memory 
effect of the spin chain reservoir. In the stationary regime, we have — Relnp(i) = 

^""2^ t- This behavior leads to exponential decay of the off-diagonal element of 
the density matrix in the long time tail. Therefore, after tracing out the spin chain 
system, the qubit density matrix exhibits a maximal mixed state, i.e., an entangled 
state is formed. In the intermediate regime (1 < ^32 < 50) , the dephasing rate 
exibits oscillatory behavior. Due to the restriction of the spin magnitude of the spin 
chain, during the transient regime between decoupled state and an entangled state 
there exists oscillatory behavior. In the stationary limit (50 -C ^r), the maximally 
entangled state is formed. Therefore, after we trace out the spin chain system, shows 
the qubit density matrix is in the maximally mixed state. The time traces qubit and 
an spin chain are separable state to maximally entangled state . 

The disappearance of the diagonal elements of the qubit indicate this behavior. 
In another context this behavior represents an effect of the qubit system on the 
environment. 

We now examine the finite temperature behavior. In the frequency regime, we 
carry out the Fourier transform of the integral in Eq. (|12|) . Then, the expression for 
the of dephasing rate is 

-Relnp(t) = ^|^ £^) { ^(0,u,),5 2 (0,-.)}((^^) 2 , (12) 
where T = 0. 

At finite temperature, because generally a spin-1/2 system obeys Fermi statistics 
we can obtain the dephasing rate by multipling the to spectral weight function by 
tanh(w/2&BT), 

-Relnp(i) = J du tanh(^-^)({5 z (0, u), S z {0, -a/)}) I — — - — 1 . 

(13) 
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The validity of the expression is demonstrated in by Ref. 17), and it should be 
required because of antiperiodic conditions (Fermi statistics). This expression is 
reminiscent of the Ingold-Nazarov theory in quantum tunneling events. 19 )~ 21 ) They 

showed the relation P(t) = expK(t), with K(t) = ± J °° ^(coth (^r) [(cos(wt)- 

1)) — i sin ujt]), where J {to) is the bath spectral density and P is a coherence func- 
tion, which represents the Coulomb-blocked effect. From this correspondence, we 
can conclude that the bath spectral density function of the Haldane-Shastry model 
at finite temperature, is given by 

J(u)hs = Yt a nHu/2k B T)\A zzlN h) 2 ({S z (0,u),S z (0,-u)}). (14) 

At finite temperatures, the spin bath has a smaller effect on the system, because of 
the possibility for saturation of the populations in the bath. A conspicuous coun- 
terexample to a spin bath is the case in which the ladder of excited states of a 
nonlinear bath mode narrows with increasing energy. In this case, we find enhance- 
ment of J{ijj) in comparison with the spectral density of a harmonic bath. This type 
of temperature dependence results from the statistics of the bath. 13 )' 15 )' 16 ) At high 
temperature spins of electrons behave like those of classical particles. However, in 
spite of this fact, such classical behavior had not been studied to this time. 

Next, we examine the situation in which dephasing with dissipation occurs 
(A± 7^ 0). In the present case, the second-order perturbation approximation is 
not accurate, although for a pure dephasing event we obtain an exact result. In 
previous studies, Fermi's Golden Rule was used. In this case, the transverse re- 
laxation, T^ 1 and the longitudinal transition rate are given by 22 ) — ln(ry = l,t) = 

Tf 1 = ^rf(3 + ^ (|) 2 ) and -ln(C = l,t) = + U A \\^K) 2 . The results 
show that the dephasing suppressed as the strength of the intrachain interaction 
increases. This is because the fluctuations of a spinon are suppressed as J increases. 
Also, as A increases, the dephasing increases, because the spin bath behaves like as 
ohmic bath. It is conjectured that for a spin bath, the higher-order corrections of 
A are negative. This behavior leads us to hypothesize a finite band width of In 
the present model, we can confirm our results, because the universality class of the 
Haldane-Shatry model is the Tomonaga-Luttinger liquid with K = 1/2, 24 ) and at 
zero temperature, the dephasing rate of dephasing with dissipation is constant. 23 ) 

§5. Summary 

In summary, we had examined the dephasing rate of a spin qubit system coupled 
with a spin chain. Due to the fluctuations of the spin of the spin chain, the dephasing 
is supperessed as J. The time dependence exhibits oscillatory behavior, which comes 
from the entangled state. We also examined dephasing with dissipation. We found 
that the dephasing with dissipation is suppressed as the strength of the magnetic 
field applied to qubit system (energy width of two-level system) is decreased. It 
should be noted that for pure dephasing, our results are exact and can be applied to 
other bath systems by using the bath spectrum density function. 
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